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Summary. We have explored in three dimensions the fate of a white dwarf of mass 
of 1.38Mq as a function of different initial locations of carbon ignition, with the 
aid of a SPH code. The calculated models cover a variety of possibilities ranging 
from the simultaneous ignition of the central volume of the star to the off-center 
ignition in multiple scattered spots. In the former case, the possibility of a tran- 
sition to a detonation when the mean density of the nuclear flame decreases to 
p ~ 2 10 7 g.cm -3 and its consequences are discussed. In the last case, the de- 
pendence of the results as a function of the number of initial igniting spots and 
the chance of some of these models to evolve to the pulsating delayed detonation 
scenario are also outlined. 



1 Statement of the Problem 

The analysis of light curves of many Type la Supernovae (SNIa) indicates that these 
events are not as homogeneous as one would desire to use them as perfect standard 
candles. In fact, the inferred mass of Nickel ejected in the explosion, estimated from 
the bolometric light curves [1], could range within a factor of ten. From a theoretical 
point of view, there remain a number of fundamental issues to be solved before we 
can rely on the predictions of SNIa modeling: pre-supernova evolution (progenitor 
systems, path of the white dwarf up to the Chandrasekhar mass, ignition condi- 
tions), physics of the explosion (subgrid-scale physics vs. large-scale combustion, 
flame behavior at low densities, deflagration-detonation transition), role of rotation 
and magnetic fields, etc. 

Although one-dimensional (spherical) models of SNIa have succeeded to explain 
a wide range of observational properties, these kind of models unavoidably rely 
on phenomenological descriptions of inherently multidimensional processes. Among 
them, flame acceleration due to hydrodynamic instabilities is of particular impor- 
tance. Recent multi(3D)dimensional simulations of deflagrations in Chandrasekhar 
mass white dwarfs have shown that the subsonic propagation of the burning front 
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can release enough energy to unbind the whole star and, maybe, can produce SNIa- 
like explosions [2, 3]. Although the final test of these models will have to be done 
by comparing self-consistent calculations of light curves and spectra (computed in 
three dimensions) with observations, one can already wonder if there is any signifi- 
cant observational evidence of departure from spherical symmetry in SNIa. In this 
respect, we note that: 1) Light curves are well described by a one-parameter family 
of curves (e.g. [4]), 2) spectral absorption features by Sill show quite homogeneous 
profiles from event to event [5], and 3) polarization has not been detected in most 
SNIa, although there are a few exceptions [6]. Supernova remnants (SNR) provide 
another means to constrain the geometry of the explosions, although it depends also 
on the presence of interstellar medium inhomogeneities and on the development of 
further hydrodynamic instabilities. Still, some SNRs do not show large departures 
from spherical symmetry (e.g. the blast wave of Tycho's SNR), which implies that 
both the supernova ejecta and the interstellar medium possessed a high degree of 
symmetry in these cases (see also the paper by Badenes et al in these proceedings). 
All these data point to approximately spherical explosions in which the chemical 
inhomogeneities are constrained to small-size clumps (a quantitative criterion, pro- 
posed in [5], is that the area of the individual clumps present at the photosphere 
has to be lower than 1 — 10% of the photospheric area). 
Our aim here is to address the following questions: 

• Are 3D simulations of thermonuclear supernovae able to cover the full range of 
light curve variations (up to a lOx factor in 56 Ni production)? 

• How much the ejecta structure obtained in these simulations deviates from spher- 
ical symmetry (overall shape, small-scale clumping)? 

We present the results from 3D hydrodynamic simulations starting from different 
initial conditions. We have explored as well the possible outcomes from deflagration- 
detonation transitions in 3D. It has to be stressed that all calculations were carried 
out consistently, i.e. using the same hydrocode [7], with the same numerical resolu- 
tion, the same physics, and simulating the whole white dwarf volume to avoid the 
introduction of artificial and unrealistic symmetry conditions. 

2 Computed 3D models, 56 Ni productivity and 
asymmetry 

We have started from several ignition configurations, in which the carbon runaway 
took place either in a central volume, slightly perturbed from spherical symmetry, 
or in multiple scattered spots or "bubbles" [8]. In the second case, we have explored 
the sensitivity of the results to questions like: How many bubbles ignite simultane- 
ously? Which are the sizes of the hot bubbles? Both theoretical analysis [9] and 2D 
simulations [10] do not seem to favor the formation of a large number of bubbles 
igniting in phase (i.e. of the same size). 

As for the deflagration-detonation transition, up to now there is no convinc- 
ing mechanism that could account for the acceleration of the burning front up to 
supersonic velocities, so we have computed delayed detonation models making use 
of several different algorithms for the initiation of the detonation. In one case, we 
started the detonation in those points in which the fractal dimension of the flame 
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Table 1. List of 3D simulations performed 



Model category Ignition configuration S kin M( 56 Ni) M(C + O) 

(10 51 erg) (M ) (M Q ) 



Deflagration 


Central volume 





.26 


0. 


.27 





.65 


Deflagration 


6 equal size bubbles 









.21 






Deflagration 


7 equal size bubbles 






0. 


.21 






Deflagration 


10 equal size bubbles 





.05 





.24 





.88 


Deflagration 


30 equal size bubbles 





.44 


0. 


.56 





.58 


Deflagration 


90 random size bubbles 





.45 


0. 


.58 





.58 


Delayed detonation 


Transition where D > 2.5 





.75 





.54 





.39 


Delayed detonation 


Transition in central layers 





.48 


0. 


.40 





.55 


Delayed detonation 


Transition in medium layers 





.51 





.43 





.51 


Delayed detonation 


Transition in external layers 





.33 


0. 


.32 





.65 


Pulsating 
















reverse detonation 


6 bubbles 





.89 





.35 





.22 



was larger than 2.5, while in the other cases, the condition for detonation ignition 
was that the density were lower than 2 x 10 7 g cm~ 3 and the particles to be det- 
onated were selected according to their radius. The initial deflagration phase of all 
the delayed detonation models was coincident with the deflagration model given in 
the first row of Table 1. 

In Table 1, there are listed the models computed and the main results of the 
simulations. In all cases, the ignition started when the white dwarf central density 
reached 1.8 x 10 9 g cm" 3 , and the simulations used 250,000 identical mass particles 
(each particle was about the same mass as the Earth). The third column in Table 1 
gives the final kinetic energy, the fourth column gives the ejected mass of 56 Ni, and 
the fifth column gives the ejected mass of unburned C and O. 

As can be seen from Table 1, all the deflagration models gave a kinetic energy 
too low to account for SNIa properties, with the exception of the two simulations 
starting from a large number (> 30) of bubbles. It is a remarkable result that the 
explosion properties resulting from these two models were almost the same, even 
though they started from very different initial conditions (in terms of number and 
sizes of hot bubbles). The results of these two simulations were as well close to the 
ones given in [3] , for comparable initial conditions, which reinforces the reliability of 
the outcome we obtained. However, in the deflagration scenario there always remains 
a large mass of unburned C and O, which would undoubtedly have to be detectable 
in the optical spectra of SNIa. 

All our delayed detonation models gave healthy explosions with reasonable ki- 
netic energies and a wide range of 56 Ni masses. If nature provided a range of condi- 
tions for transition to detonation as wide as we have explored, the delayed detonation 
scenario would help to explain the diversity of SNIa explosions observed. However, 
our 3D simulations have shown that the delayed detonations are relatively inefficient 
in converting C and O to intermediate mass elements or to iron peak elements (which 
is in sharp contrast with ID models of delayed detonations), which results in a lower 
velocity range than expected. The ultimate reason for this behavior is the geometry 
of the distribution of fuel resulting from the previous deflagration stage. Due to the 
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Fig. 1. Final distribution of unburned C-0 in selected models. From left to right: 
deflagration starting from a central volume, deflagration starting from 30 bubbles, 
and delayed detonation initiated where the flame fractal dimension was D > 2.5 



unavoidable distortion of the flame front caused by hydrodynamic instabilities dur- 
ing subsonic burning, the initial conditions for delayed detonation formation consist 
on large plumes of ashes coexisting at the same radii with deep tongues of fuel (see 
Fig. 1). These ashes' plumes act as true barriers obstructing the propagation of the 
detonation waves and, in fact, decoupling large volumes occupied by C and O from 
other regions which can hold propagating detonations. 

To discuss the geometrical properties of the explosions, we will refer to Fig. 
1, in which we show the distributions of the ejected C-O in YZ slices. None of 
the models displays an overall shape that departs largely from spherical symmetry, 
although the distribution of chemical species presents irregularities of different sizes. 
However, the deflagration model starting from a central volume (left image in Fig. 
f ) is the most asymmetric. This is due to Rayleigh-Taylor instabilities during flame 
propagation, which favor the formation of large-scale structures that stand out in 
the final geometry of the ejecta. This configuration is in clear disagreement with the 
limits to the size of clumps given in [5] . In the case of deflagrations starting from hot 
bubbles, the geometrical appearance of the ejecta retains a larger degree of spherical 
symmetry, but still the chemical inhomogeneities represent a large fraction of the 
radius of the star. We have estimated that, in the model starting from 30 bubbles 
(center image in Fig. f), a typical clump size is ~ 10% of the ejecta radius, still too 
large. 

Delayed detonations (right image in Fig. 1) are approximately spherical in shape, 
and show only small clumps. We have estimated a typical clump size of ~ 4% 
of the ejecta radius in these kind of models. This could seem a striking result, 
as the delayed detonations were computed starting from the final configuration of 
the deflagration shown in the left image in Fig. 1. The reason for this apparent 
contradiction is that the detonations were able to destroy the large structures built 
up during the deflagration phase, thus allowing the formation of only small-sized 
chemically differentiated clumps. 



3 Discussion and a new SNIa paradigm 

To summarize, deflagrations computed in 3D produce large clumps, ejecta with 
low kinetic energy, reasonable quantities of 56 Ni, and too much unburned C-0 (> 
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0.58 M Q ). On the other hand, multidimensional delayed detonations produce smaller 
clumps, although they are not as efficient as their ID counterparts at rising the 
kinetic energy (< 0.75 x 10 51 erg) nor at burning C and O (M{C - O) > 0.39 M©). 
In addition, we have seen that deflagrations starting from a reasonable number of 
bubbles (< 6 — 10) fail to unbind the white dwarf, which brings us to the following 
question: What is the fate of the white dwarf following explosion failure? 

We have followed the evolution of the star during the first pulsation with the 
3D SPH code. What came out from this multidimensional calculation was in fact 
quite different from what is obtained in ID pulsating models. Due to the ability of 
bubbles to float to large radii in 3D, most of the thermal and kinetic energy resided 
in the outer parts of the structure, which resulted in an early stabilization of the 
central region (mostly made of C and O, i.e. fuel) while the outer layers were still in 
expansion. Few seconds later, an accretion shock formed at the border of the central 
nearly hydrostatic core (whose mass was about 0.9 M ). Therefore, the temperature 
at the border of the core increased to nearly 10 9 K, on a material composed mainly 
by fuel but with a non-negligible amount of hot ashes, thus giving rise to a highly 
explosive scenario. If a detonation were ignited at this point, it would probably 
propagate all the way inwards through the core, burning most of it and producing an 
energetic explosion. We have called this new paradigm of SNIa explosion mechanism 
the Pulsating Reverse Detonation. Currently, we are in the process of completing 
the 3D calculation of its detonating phase. As a first evaluation of what one can 
expect from this new scenario, we have computed a ID model (last row in Table 1), 
which ended with the largest kinetic energy, and the lowest mass of C and O from 
all the set of models we have computed in 3D up to now. It turns out that it is 
worth following the evolution of this scenario in 3D to see if it can derive in a new 
reliable and competitive model of SNIa. 
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